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Abstract. We apply the finite size scaling technique to 
quantify the statistical properties of fluctuations in AU, AL 
and AE indices and in the e parameter that represents en- 
ergy input from the solar wind into the magnetosphere. We 
find that the exponents needed to rescale the probability 
density functions (PDF) of the fluctuations are the same to 
within experimental error for all four quantities. This self- 
similarity persists for time scales up to ~ 4 hours for AU, AL 
and e and up to ~ 2 hours for AE. Fluctuations on shorter 
time scales than these are found to have similar long-tailed 
(leptokurtic) PDF, consistent with an underlying turbulent 
process. These quantitative and model-independent results 
place important constraints on models for the coupled solar 
wind-magnetosphere system. 



1. Introduction 

Recently, there has been considerable interest in viewing 
the coupled solar wind-magnetosphere as a complex system 
where multi-scale coupling is a fundamental aspect of the 
dynamics (see [Chang, 1992; Chapman and Watkins, 2001] 
and references therein). Examples of the observational mo- 
tivation for this approach are i) bursty transport events in 
the magnetotail [Angelopoulos et at, 1992] and ii) evidence 
that the statistics of these events are self-similar (as seen 
in auroral images [Lui et al, 2000]). Geomagnetic indices 
are of particular interest in this context as they provide a 
global measure of magnetospheric output and are evenly 
sampled over a long time interval. There is a wealth of 
literature on the magnetosphere as an input-output system 
(see for example, [Klimas et al, 1996; Sitnov et at, 2000; 
Tsurutani et al., 1990; Vassiliadis et al., 2000; Voros et al, 
1998]). Recent work has focussed on comparing some as- 
pects of the scaling properties of input parameters such as 
e [Perreault and Akasofu, 1978] and the AE index [Davis 
and Sugiura, 1966] to establish whether, to the lowest or- 
der, they are directly related [Freeman et al, 2000; Uritsky 
et al., 2001]. Although these studies are directed at un- 
derstanding the coupled solar wind-magnetosphere in the 
context of Self-Organized Criticality (SOC), a comprehen- 
sive comparison of the scaling properties of the indices, and 
some proxy for the driver (e) also has relevance for the pre- 
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dictability of this magnetospheric "output" from the input. 
Importantly, both "burstiness" (or intermittency) and self- 
similarity can arise from several processes including SOC 
and turbulence. Indeed, SOC models exhibit threshold in- 
stabilities, bursty flow events and statistical features consis- 
tent with the "scale-free" dynamics such as power law power 
spectra. It has been proposed by Chang [1992] that magne- 
tospheric dynamics are indeed in the critical state or near 
it. Alternatively, Consolini and De Michelis [1998] used the 
Castaing distribution - the empirical model derived in Cas- 
taing et al. [1990] and based on a turbulent energy cascade 
- to obtain a two parameter functional form for the Prob- 
ability Density Functions (PDF) of the AE fluctuations on 
various temporal scales. Turbulent descriptions of magne- 
tospheric measures also model observed statistical intermit- 
tency, i.e., the presence of large deviations from the aver- 
age value on different scales [Consolini et al, 1996; Voros 
et al, 1998]. An increased probability of finding such large 
deviations is manifested in the departure of the PDF from 
Gaussian toward a leptokurtic distribution [Sornette, 2000]. 

In this paper we will quantify both the intermittency and 
the self-similarity of the AU, AL, AE and e time series using 
the technique of finite size scaling. This has the advantage 
of being model independent, and is also directly related to 
both turbulence models such as that of Castaing [Castaing 
et al, 1990] and a Fokker-Planck description of the time se- 
ries. The method was used in Hnat et al. [2002] where the 
mono-scaling of the solar wind magnetic energy density fluc- 
tuations was reported. We will find that fluctuations in all 
four quantities are strongly suggestive of turbulent processes 
and by quantifying this we can compare their properties di- 
rectly. 

The AL, AU and AE indices data set investigated here 
comprises over 0.5 million, 1 minute averaged samples from 
January 1978 to January 1979 inclusive. The e parameter 
defined in SI units as: 



where lo ~ 7 'Re and = axctan(\B v \/B z ) is an estimate 
of the fraction of the solar wind Poynting flux through the 
dayside magnetosphere and was calculated from the WIND 
spacecraft key parameter database [Lepping et al., 1995; 
Ogilvie et al., 1995]. It comprises over 1 million, 46 sec- 
ond averaged samples from January 1995 to December 1998 
inclusive. The data set includes intervals of both slow and 
fast speed streams. The time series of indices and that of 
the e parameter were obtained in different time intervals 
and here we assume that the samples are long enough to be 
statistically accurate. 

2. Scaling of the indices and e 

The statistical properties of complex systems can exhibit 
a degree of universality reflecting the lack of a characteristic 
scale in their dynamics. A connection between the statistical 
approach and the dynamical one is given by a Fokker-Planck 
(F-P) equation [van Kampen, 1992] which describes the dy- 
namics of the PDF and, in the most general form, can be 
written as: 



e = v — llsm 4 (Q/2), 



(1) 



dP(x,t) 
dt 



= V(P(x,t) 1 (x)) + V 2 D(x)P(x,t), (2) 



!! Please write \lef thead{<AUTHOR NAME(s)>} in file !!: !! Please write \righthead{< (Shortened) Article Title>} in file !! 3 



where P(x, t) is a PDF of some quantity x that varies with 
time t, 7 is the friction coefficient and D(x) is a diffusion 
coefficient which in this case can vary with x. For certain 
choices of D(x), a class of self-similar solutions of (2) satisfies 
a finite size scaling (in the usage of Sornette [2000], pg. 85, 
henceforth "scaling") relation given by: 

P(x,t) =t- s P a (xt- s ). (3) 

This scaling is a direct consequence of the fact that the F-P 
equation is invariant under the transformation x — ► xt~ s and 
P — ► Pt 3 . ff, for given experimental data, a set of PDFs can 
be constructed, on different temporal scales r, that satisfy 
relation (3) then a diffusion coefficient and corresponding 
F-P equation can be found to represent the data. A sim- 
ple example is the Brownian random walk with s = 1/2, 
D(x) ^constant and Gaussian PDFs on all scales. Alter- 
natively one can treat the identification of the scaling ex- 
ponent s and, as we will see, the non-Gaussian nature of 
the rescaled PDFs (P s ) as a method for quantifying the in- 
termittent character of the time series. Practically, obtain- 
ing the rescaled PDFs involves finding a rescaling index s 
directly from the integrated time series of the quantity X 
[Hnat et ai, 2002; Sornette, 2000]. 

Let X(t) represent the time series of the studied signal, 
in our case AU, AL, AE or the e parameter. A set of time 
series SX(t, r) = X(t + r) — X(t) is obtained for each value 
of non-overlapping time lag r. The PDF P(SX, r) is then 
obtained for each time series 5X(t, r). Figure 1 shows these 
PDFs for the SAU. A generic scaling approach is applied to 
these PDFs. Ideally, we use the peaks of the PDFs to obtain 
the scaling exponent s, as the peaks are the most populated 
parts of the distributions. In certain cases, however, the 
peaks may not be the optimal statistical measure for ob- 
taining the scaling index. For example, the B z component 
in (1) as well as the AU and AL indices are measured with 
an absolute accuracy of about 0.1 nT. Such discreteness in 
the time series and, in the case of the e fluctuations, the 
large dynamical range introduce large errors in the estima- 
tion of the peak values P(0, r) and may not give a correct 
scaling. Since, if the PDFs rescale, we can obtain the scaling 
exponent from any point on the curve in principle, we also 
determine the scaling properties of the standard deviation 
<j(t) of each curve P(SX,t) versus SX(t,r). 

Figure 2 shows P(0, r) plotted versus r on log-log axes 
for SX = Se, SAE, SAU and SAL. Straight lines on such 
a plot suggest that the rescaling (3) holds at least for the 
peaks of the distributions. On figure 2, lines were fitted with 
R 2 goodness of fit for the range of r between 4 and 136 min- 
utes, omitting points corresponding to the first two temporal 
scales as in these cases the sharp peaks of the PDFs can not 
be well resolved. The lines suggest self-similarity persists up 
to intervals of r = 97 — 136 minutes. The slopes of these 
lines yield the exponents s and these are summarized in Ta- 
ble 1 along with the values obtained from analogous plots 
of ct(t) versus r which show the same scale break. We note 
that, for the e parameter, the scaling index s obtained from 
the P(0, t) is different from the Hurst exponent measured 
from the cr(r). This difference could be a result of the previ- 
ously discussed difficulties with the e data. However, it does 
appear to be a feature of some real time series (see Gopikr- 
ishnan et al [1999] for example). Indeed, such a difference 
between index s and H a is predicted in the case of fractional 
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Levy motion [Chechkin and Gonchar, 2000]. We see that, 
for the e as well as the AL and AU indices, there is a range 
of t up to 4.5 hours for which P(0, r) plotted versus r is well 
described by a power law r~ s with indices s = 0.42 ± 0.03 
for the e and s = 0.45 ± 0.02 and s = 0.47 ± 0.03 for the AL 
and AU indices, respectively. Thus the break in scaling at 
4 — 5 hours in the AL and AU indices may have its origin in 
the solar wind, although the physical reason for the break 
at this timescale in epsilon is unclear. The break in the AE 
index, however, appears to occur at a smaller temporal scale 
of 2 hours, consistent with the scaling break timescale found 
in the same index by other analysis methods [Consolini and 
De Michelis, 1998; Takalo et al, 1993]. This was interpreted 
by [Takalo et al, 1993] as due to the characteristic substorm 
duration. Takalo and Timonen [1998] also reported a scal- 
ing break at the same 2 hour timescale for AL, in contrast 
to the 4 — 5 hour timescale found here. Indeed, one might 
have expected a substorm timescale to cause the same scal- 
ing break in both the AE and AL indices, because their 
substorm signatures are so similar in profile (e.g., Figure 2 
of Caan et al. [1978]). The resolution may lie in the differ- 
ence between analysis of differenced and undifferenced data 
[Price and Newman, 2001]. 

Within this scaling range we now attempt to collapse each 
corresponding unsealed PDF onto a single master curve us- 
ing the scaling (3). If the initial assumption of the self- 
similar solutions is correct, a single parameter rescaling, 
given by equation (3) for a mono-fractal process, would give 
a perfect collapse of PDFs on all scales. Practically, an ap- 
proximate collapse of PDFs is an indicator of a dominant 
mono-fractal trend in the time series, i.e., this method may 
not be sensitive enough to detect multi-fractality that could 
be present only during short time intervals. Figures 3 and 
4 show the result of the one parameter rescaling applied to 
the unsealed PDF of the Se and the SAU index fluctuations, 
respectively, for the temporal scales up to ~ 4.5 hours. We 
see that the rescaling procedure (3) using the value of the 
exponent s of the peaks P(0, r) shown in Fig 2, gives good 
collapse of each curve onto a single common functional form 
for the entire range of the data. These rescaled PDFs are 
leptokurtic rather than a Gaussian and are thus strongly 
suggestive of an underlying turbulent process. 

The successful rescaling of the PDFs now allows us to 
perform a direct comparison of the PDFs for all four quan- 
tities. Figure 5 shows these normalized PDFs P s (SX,r) for 
SX = Se, SAE and r » 1 hour overlaid on a single plot. The 
SX variable has been normalized to the rescaled standard 
deviation a 3 (r ~ lhr) of P 3 in each case to facilitate this 
comparison. We then find that AE and e fluctuations have 
indistinguishable P s . The PDFs of SAU and SAL are asym- 
metric such that —SAL fits SAU PDF closely (see insert in 
the figure 5); when overlaid on the PDFs of the Se and SAE 
these are also indistinguishable within errors. This provides 
strong evidence that the dominant contributions to the AE 
indices come from the eastward and westward electrojets 
of the approximately symmetric DP2 current system that 
is driven directly by the solar wind [Freeman et al, 2000]. 
The mono-scaling of the investigated PDFs, together with 
the finite value of the samples' variance, indicates that a 
Fokkcr-Planck approach can be used to study the dynamics 
of the unsealed PDFs within their temporal scaling range. 
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3. Summary 

In this paper we have applied the generic and model inde- 
pendent scaling method to study the scaling of fluctuations 
in the e parameter and the global magnetospheric indices 
AU, AL and AE. The similar values of the scaling exponent 
and the leptokurtic nature of the single PDF that, to within 
errors, describes fluctuations on time scales up to r max in all 
four quantities provide an important quantitative constraint 
for models of the coupled solar wind-magnetosphere system. 
One possibility is that, up to T ma x ~ 4 hours, fluctuations in 
AU and AL are directly reflecting those seen in the turbu- 
lent solar wind. The data also suggest that AE index departs 
from this scaling on shorter time scale of r max ~ 2 hours. 
Importantly, identifying a close correspondence in the fluc- 
tuation PDF of e, AE, AU and AL may simply indicate that 
fluctuations in the indices are strongly coupled to dayside 
processes and are thus weak indicators of the fluctuations 
in nightside energy output. The leptokurtic nature of the 
PDFs is strongly suggestive of turbulent processes, and in 
the case of AU and AL, these may then be either just that of 
the turbulent solar wind (and here e) or may be locally gen- 
erated turbulence which has an indistinguishable signature 
in its fluctuation PDF. In this case our results quantify the 
nature of this turbulence. We note, however, that certain 
classes of complex systems [Chang et ai, 1992] are in prin- 
ciple capable of "passing through" input fluctuations into 
their output without being directly driven in the present 
sense [Chang, private communication, 2002]. Finally, the 
rescaling also indicates that a Fokker-Planck approach can 
be used to study the evolution of the fluctuation PDF. This 
raises a possibility of a new approach to understanding mag- 
netospheric dynamics. 
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Figure 1. Unsealed PDFs of the AU index fluctuations. 
Time lag r assumes values between 60 seconds and about 36 
hrs. Standard deviation of the PDF increases with r. Error 
bars on each bin within the PDF are estimated assuming 
Gaussian statistics for the data within each bin. 
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Figure 2. Scaling of the peaks of the PDFs for all quantities 
under, investigation: o corresponds to e, o AU index, A 
AL index and □ the AE index. The plots have been offset 
vertically for clarity. Error bars as in Figure 1. 
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Quantity 


P(0, t) scaling 
index 


ct(t) scaling 
index 


Tvaax 


e 


-0.42 ± 0.03 


0.33 ± 0.04 


4.5 hrs 


AE-index 


-0.44 ±0.03 


0.43 ± 0.03 


2.1 hrs 


AU-indcx 


-0.47 ± 0.03 


0.47 ±0.02 


4.5 hrs 


AL-index 


-0.45 ± 0.02 


0.45 ± 0.02 


4.5 hrs 



Table 1. Scaling indices derived from P(0,r) and <t(t) power 
laws. 
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Figure 3. One parameter rescaling of the e parameter fluc- 
tuations PDFs. The curves shown correspond to r between 
46 seconds and 4.5 hours. Error bars as in Figure 1. 
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Figure 4. One parameter rescaling of the AU index fluctu- 
ation PDF. The curves shown correspond to r between 46 
seconds and 4.5 hours. Error bars as in Figure 1. 
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Figure 5. Direct comparison of the fluctuations PDFs for e 
(o) and AE index (□). Insert shows overlaid PDFs of AU(o) 
and — AL(A) fluctuations. Error bars as in Figure 1. 



